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ABSTR/LCT 


A basic  analysis  of  the  blade  tip  flow  phenomena  is  presented# 
The  analysis  is  based  on  a steady  relative  flow  of  a viscous  incompres- 
sible fluid.  A rotating  cylindrical  coordinate  system  which  moves 
with  the  blade  at  the  same  speed  is  used. 

An  order-of-magnitude  analysis  reduces  the  basic  equations  to 
a set  of  simplified  equations  for  both  low  and  high  Reynolds  number# 
Expressions  are  given  for  the  determination  of  the  distributions  of 
pressure  and  velocity  in  the  clearance  space  and  the  mass  flow  across 
the  tip  clearance  space.  The  computation  of  the  velocity  distribution 
is  relatively  single  and  three  methods  of  solution  for  the  pressure 
distribution  are  presented# 
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INTRODUCTION 

Because  of  the  relative  potion  between  the  rotating  blades 

and  the  stationary  casing  wall  in  the  turbomachines,  a clearance  space 

has  to  be  provided  to  avoid  physical  rubbing  and  consequent  damaging 

of  the  machine.  The  existence  of  the  tip-clearance  complicates 

greatly  the  internal  fluid  flow  in  turbomachines.  Briefly,  the  tip- 

clearance  flow  is  a three-dimensional  flow  problem.  The  tip-clc franco 

flow  has  long  been  considered  to  be  one  of  the  factors  responsible  for 

(l)* 

internal  losses  in  trubomachines,  although  opinion  of  its  importance 
varies.  The  loss  due  to  tip  clearance  flow  has  also  been  known  to  be 
closely  associated  with  the  secondary-flow  loss  and  the  annulus-wall 
loss.  The  tip-leakage  flow  is  due  to  the  pressure  difference  across 
the  blade-tip  section  and  the  relative  motion  between  the  blade  and 
the  bounding  wall  and  is  transverse  to  the  main-flow  direction.  The 
secondary-flow  is  due  to  the  pressure  difference  across  the  space 
between  two  adjacent  blades  and  is  also  transverse  to  the  main-flow 
direction.  The  annulus  wall  loss  is  due  to  the  boundary  layer  at  the 
wall.  All  these  losses  are  originated  within  the  boundary  layers  of 
the  casing  wall  surfaces. 


Of  course,  the  tip-clearance  should  be  kept  as  small  as  possible 

But  the  "warm  clearance"  varies  over  a wide  range  under  different 

(2) 

operating  conditions.  It  is  therefore  advisable  to  provide  a liberal 
clearance  in  order  to  avoid  rubbing  under  most  unfavorable  conditions, 
"Shrouding"  of  blade  tips  decreases  the  tip-loss  but  introduces  other 

dumber  in  parenthesis  refers  to  the  reference  number  given  at  the  end 
of  the  report. 
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(3)  (W 

losses.  In  one  compressor  test,  it  has  been  found  that  the  compressor 
with  shrouded  blades  has  both  lower  pressure  coefficient  and  overall 
compressor  efficiency  than  those  of  a compressor  without  shrouding. 

(5) 

However,  this  single  test  cannot  be  taken  as  conclusive.-  Others 
claim  shrouding  will  eliminate  the  induced  secondary  flow  loss  which  is 
caused  by  the  moving-wall.  It  seems  that  shrouding  only  changes  the 
tip-loss  problem  to  a different  one  instead  of  being  a simple  remedy. 
Our  study  is  confined  to  free  tip  blades. 

The  purpose  of  this  investigation  is  to  make  a basic  study  of 
the  flow  phenomenon  pattern  around  the  tip-clearance  space,  to  analyze 
the  design  and  operating  parameters  which  control  this  clearance  flow, 
and  to  determine  its  effect  on  the  machine  performance.  It  is  hoped 
that  out  of  a thorough  understanding  of  the  tip-flow  phenomena,  the 
limit  of  the  permissible  clearance  may  be  calculated.  This  report 
presents  the  initial  theoretical  study  of  this  problem. 

A BRIEF  SURVEY  OF  PREVIOUS  INVESTIGATIONS 

Up  to  date,  complete  treatment  of  the  tip-loss  problem  is  not 

available,  even  though  it  has  attracted  attention  in  the  steam-turbine 

(6) 

development  as  early  as  190$,  The  available  information  regarding 
the  tip-loss  problem  can  be  summarized  into  the  following  different 
groups: 

Theoretical  Investigation 

(7) 

Betz  made  the  first  theoretical  investigation  of  the  tip-loss 

in  the  Kaplan  water  turbines  by  applying  the  lifting  line  theory  to  a 

(8) 

simplified  two-dimensional  rectilinear  cascade.  Later,  Sedille 
also  using  lifting-line  theory,  treated  the  tip-loss  problem  in  the 
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axial  compressor.  The  results  obtained  by  this  method  were  found 
much  larger  than  those  observed  in  practice.  The  discrepancy  was 
attributed  by  Sedille  to  the  resistances  which  opposed  the  flow  in 
the  tip  clearance  space  which  prevails  in  the  actual  machines  but 
are  neglected  in  the  analyses.  These  resistances  included  annulus 
wall  boundary  layers,  blade  thickness,  radial  flow  of  blade  boundary 
layers,  local  turbulence,  etc. 

Semi-Empirical  Estimation  of  Tip-Loss 


In  order  to  meet  the  urgent  practical  need  in  design,  various 

semi-cmpirical  formulae  have  been  suggested  to  estimate  the  tip-loss 

in  turbomachines.  These  formulae  were  chiefly  based  upon  simplified 

theoretical  analyses  under  certain  assumptions  or  limitations.  The 

empirical  constants  have  to  be  evaluated  from  tost  data.  This  group 

(9)  (10)  (11) 

of  investigators  includes  Meldahl,  Traupel  and  Fickert.  Their 
suggested  formulae  respectively  are; 


Investigators 

(1)  Meldahl (19^1) 

VJhere 


Suggested  Formulae 
End  Losses: 

Lq  - .1011  + U.667(-5-) 

c 

s ■ tip-clearance 
c * blade -chord 


Assumptions  and 
Limitation 


(1)  Based  upon 
lift-line  theory 
and  2 -Dimensional 
blades  cascade 

(2)  Good  only  for 
single  stage 
reaction  turbine 

(3)  Including 
secondary  flow 
loss 
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(2)  Traupel  (19U2) 


Tip  Loss:  > 

(1) 

Eased  on  Bernoulli* s 

\ 

equation 

Lt  - K(Ac 

(2) 

Flow  coeff.  * K*  has  to 

Aa  sin  o' 

be  evaluated  from  test 

data 

Where :K  * Flow  Coeff,  (3)  K has  different  values 

for  rotor  and  stator 

Ac  “ clearance  area 
Aa  “ annular  area 
t ■ pitch 

AZJ  ■ change  of  tangen- 
tial component  of 
velocity 

■ cascade  angle 


(3)  Fickert  (19U3)  Volumetric  eff. 


(1)  Based  upon  Bernoulli* s 
equation 


7 


Coeff  of  contraction 
*/( * has  to  be  evalua- 
ted from  test  data 


Where: n*  coeff.  of  contrac- 
/ tion 

Da  **  casing  diameter 
= hub  " 

<T  » throttling  coeff. 


Experimental  Investigation 


In  axial  compressor  tests  the  variation  of  compressor  efficiency 

(8)  (12)  (13) 

with  tip-clearance  was  obtained  by  Sedillc,  Ruden  and  Lindsay, 

These  results  show  a nearly  linear  relationship  betweon  the  efficiency 
drop  and  the  tip-clearance  increase,  Howovor,  neither  detailed 
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explanations  nor  theory  were  advanced  from  those  testing  results,* 

Qualitative  Investigation  of  the  Moving-Wall  Effect 
Recently,  efforts  have  been  made  to  understand  the  moving-wall 
effect  on  the  clearance  flow.  This  is  an  at; tempt  to  separate  one  fac- 

(5) 

tor  from  others  affecting  tho  clearance  flow.  Ainley  and  Jeffs 
concluded  from  their  compressor  tost  that  the  rotating  drum  drags  along 
the  adjacent  fluid  through  the  tip-clearance  space  in  addition  inducing 
a secondary  flow.  Carter  indicated  that,  in  the  case  of  a compres- 
sor, tho  clearance  flow  will  bo  augmented  by  the  moving  wall  and  that 
the  induced  secondary  flow  was  the  "scraping  effect ‘ of  the  moving-wall; 
but  in  the  turbine,  the  clearance  flow  will  be  reduced  by  the  moving 
wall.  Both  the  * scraping  effect1  and  the  reduced  clearance  flow  in 

(15) 

the  turbine  wore  confirmed  by  Hansen,  Hcrzig,  Costello,  in  their  smoke 
visualization  oxperimont  in  a cascade  tunnel.  But  their  results 
regarding  the  clearance  flow  in  the  compressor  were  in  contradiction 
to  Carter* s prediction.  Their  experiments  show  that  the  moving  wall 
has  a tendency  to  diminish  the  clearance  flow  rather  than  to  promote  it. 
This  point  needs  clarification  by  further  experiments,  and  will  be  one 
of  the  objects  to  bo  studied  in  the  experimental  part  of  this  research 
project. 


In  hydraulic  turbomachincs,  such  as  pumps  and  marine  ducted 

(16) 

propellers,  the  so-called  "vortox  cavitation"  has  been  observed. 

Fluid,  passing  through  the  tip  clearance,  forms  the  tip  vortex.  At  tho 

#it  has  come  to  our  attention  very  recently  that  a report  Tip-Loss 
has  been  published  in  Germany  in  1952,  However,  that  report  is  not 
yet  available  to  us  at  the  present  time. 
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center  of  this  vortex,  cavities  arc  produced.  These  cavities  appear 
and  form  a nearly  continuous  core  which  starts  at  about  tho  leading 
edge  of  tho  blado  and  extends  downstream.  The  magnitude  of  the 
cavities  seems  to  depend  upon  the  tip-clearance. 

From  the  preceding  survey  of  the  previous  investigations,  wo  may 
conclude  that: 

(1)  Tip-cloaranco  flow  contributes  considerable  loss  in  turbo- 
machines, especially  when  the  clearance  is  largo, 

(2)  Tip-cloaranco  flow  is  affected  by  factors,  such  as  blade- 
tip  thickness,  blade  loading  (especially  loading  conditions  near  tip 
section),  relative  motion  effect,  annulus  wall  boundary  layer,  blado 
boundary  layer,  etc.  Therefore,  it  is  a 3-D*al  boundary  layer  problem. 

PRELIMINARY  STUDY  OF  TIP  LOSS  PROBLEM 

When  a lifting  surface  (either  a wing  or  a blade)  is  placed  in 

line  with  a passing  fluid,  a true  2-D’al  flow  can  be  achieved  only  if 

the  lifting  surface  has  an  infinite  span.  However,  in  reality,  the 

lifting  surface  has  only  finite  span,  therefore,  3-D‘al  flow  phenomena 

prevails.  In  the  case  of  a wing,  the  lifting  surface  has  a free  end,  or 

tip.  The  fluid  on  the  pressure  surface  tries  to  flow  to  the  suction 

surface  over  the  free  tip.  Consequently,  a trailing  vortex  sheet  is 

formed  behind  tho  wing,  (Fig.  1),  The  vortox  sheet  causes  the  induced 

drag  of  tho  wing.  In  the  case  of  a cascado  of  blades  which  are  attached 

to  end  walls,  secondary  flow  is  developed  both  in  the  end  wall  boundary 

(17)  (18) 

layer  and  in  the  main  flow  (passage  vortex).  The  boundary  layer  fluid 

on  the  pressure  surface  of  ono  blado  flows  to  tho  suction  surface  of  its 

adjacent  blade  along  the  end  wall,  (Fig.  2).  Downstream  of  the  cascado, 
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a trailing  vortex  sheet  again  is  formed,  If  "clearance"  exists  between 
the  blade  tip  and  the  end  wall  (consider  for  the  time  being  that  the 
end  wall  is  stationary) , both  the  secondary  flow  in  the  end  wall  and 
the  trailing  vortices  will  still  be  present  but  with  different  magnitude* 
Because  of  the  clearance,  the  fluid  on  the  prossuro  surface  will  tend 

i 

to  flow  to  the  suction  surface  of  the  same  blade  (since  the  passage  is 
shorter  and  the  prossure  difference  is  the  same)  instead  of  to  the 
suction  surface  of  the  adjacent  blade.  (Fig.  3).  Due  to  this 
clearance  flow,  the  pressure  gradient  between  the  pressure  and  the 
suction  surfaces  will  no  longer  be  the  same  as  that  when  there  is  no 
clearance.  Therefore,  the  secondary  flow  will  somewhat  decrease  in 
strength.  No  doubt,  this  rearrange  of  flow  pattern  in  the  boundary 
xayer  will,  in  turn,  influence  the  flow  in  the  main  stream.  It  may  be 
noted  here  that  both  the  secondary  flow  and  the  clearance  flow  are  in 
a plane  transverse  to  the  main-flow  direction. 

Now,  let  the  end  wall  move,  so  we  can  have  a condition  corres- 
ponding to  the  actual  one  in  turbomachines.  In  the  compressor,  the 
pressure  surface  is  loading  in  the  moving  direction  and  in  the  turbine, 
the  suction  surface  is  leading.  For  the  time  being,  lot  us  agreo  that 
the  moving  wall,  duo  to  viscous  effect,  will  drag  along  the  fluid  with 
it  when  it  passes  over  the  blade-tip.  Thus,  in  the  compressor,  because 
the  wall  (rclativo  motion)  moves  in  the  same  direction  to  that  of  the 
pressure  difference  (hence  in  the  direction  of  the  clearance  flow) 
across  the  blade  tip,  it  is  expected  that  the  clearance  flow  will  be 
augmented  (Fig.  it)  while  in  the  turbine,  the  wall  moves  in  an  opposite 
direction  to  that  of  the  pressure  difference,  it  is  expected  that  tho 
clearance  flow  will  be  reduced  (Fig.  £), 
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it seems  that  it  is  necessary  to  point  out  that  there  i3  a 
difference  botwcon  the  amount  of  fluid  flow  through  the  tip-clearance 
and  the  strength  of  the  vortex  formed  due  to  this  tip-clearance  flow. 
The  amount  of  fluid -flow  through  the  tip-clearance,  as  will  be  seen 
later,  depends  upon  both  the  pressure  gradiont  across  the  blade  tip- 
section  and  the  moving  wall  velocity.  But  the  vortex  strength  depends 
only  upon  the  pressure  gradient  across  the  blade  tip-section.  Thore 
are  two  extreme  cases:  (a)  For  stationary  wall  (zero  velocity)  and 

large  prossuro  gradiont,  the  tip-clcaranco  vortex  will  be  the  strongest, 
(b)  For  largo  moving  wall  velocity  and  zero  pressure  gradient,  there 
will  be  no  tip-clearance  vortex  and  the  amount  of  fluid  flow  may  be  or 
may  not  be  greater  than  that  in  case  (a) . In  an  actual  machine,  the 
situation  for  the  vortex  strength  and  the  amount  of  fluid  flow  will  be 
somewhere  between  cases  (a)  and  (b) . The  moving  wall,  of  course,  tends 
to  minimize  tho  pressure  gradiont  across  the  two  sides  of  the  blade 
when  it  drags  the  high-pressure  fluid  from  one  side  to  the  low  pres- 
sure region  of  tho  other  side  of  the  blade. 

Because  the  flow  in  the  boundary  layer  is  primarily  determined 
by  the  main  flow,  all  the  factors  which  affect  tho  main  flow  will  affect 
the  flow  in  the  boundary  layer.  The  boundary  layer  flow  condition  in 
turn  affects  the  tip-clearance  flow.  Therefore,  tho  tip  clearance  flow 
is  closely  rolatcd  to  the  main  flow.  The  final  forms  of  the  main  flow 
influence  are  manifested^  by  the  blade  loading,  casing  wall  boundary 
layer,  relative  motion  effect,  etc.,  as  mentioned  before.  If  these 
conditions  prevailing  outside  tho  clearance  space  remain  the  same,  tho 
flow  through  the  tip-clearance  can  bo  expected  to  be  the  same  even 
though  tho  main  flow  could  be  of  different  patterns 
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for  different  designs.  Since  the  main  flow  in  the  turbomachine  has 

(19) 

many  natures  on  the  one  hand,  and  the  lack  of  full  knowledge 
regarding  the  throe  dimensional  boundary  layer  flow  (in  this  case,  it 
will  be  secondary  flow)  on  the  other  hand,  it  seems  unwise  to  attempt 
to  solve  the  tip-clearance  problem  together  with  the  boundary  layer 
flow  and  the  main  flow,  but  rather  to  assume  conditions  outside  the 
clearance  space  and  try  to  solve  it  accordingly.  In  other  words,  the 
boundary  conditions  of  tip-clearance  flow  are  not  evident  by  them- 
selves. They  are  rolated  to  the  secondary  flow  and  hence  the  main  flow. 

So  far  we  consider  the  blade  tip  is  inside  the  casing-wall 
boundary  layer.  In  other  words,  we  have  considered  only  the  case  that 
the  "warm  clearance"  is  less  than  the  thickness  of  the  boundary  layer 
of  the  adjacent  casing  wall.  This  will  be  true  for  most  conditions,  at 
least  for  later-stage  rotors  and  stators  in  the  compressor.  In  the 
following,  we  will  concentrate  our  analysis  on  this  assumption, 
Fortunately,  for  the  caso  when  the  blade  tip  is  outside  of  the  boundary 
layer,  the  problem  is  much  simpler  to  handle.  Because  both  viscous 
effect  and  the  relative  motion  effect  can  be  ignored  when  the  tip  is 
outside  of  the  boundary  layer  (Fig.  6),  then  the  problem  can  bo 
approached  by  the  conventional  lifting-line  thoory. 
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SYMBOLS 


The  following  symbols  are  used  in  this  report: 


A » Coefficient  of  cosines  of  the  velocity  Fourier  series 
B * Coefficient  of  sines  of  the  velocity  Fourier  series 
C • Airfoil  chord  length  - ft, 

D *(r3  j.  + £<r-F  7».v/3  V ) + (zef  )(&€*;*+  +£  rAN/ZMlf- 

4-  c TAM * ft)  + c*  7>tcX '($  it2'  --it 

L * Characteristic  length  of  the  turbomachine,  cither  the  blade  chord 
or  the  blade  height  - ft* 

M • l4(l+f)c  R TAN  f*  (V  )]*■'*  V-  'jty'-  fR  ( & "'-?/?  D/oz-iy 

+(-§S  ~ <-  *0  ~£<P  ~Cz)c*t  *■  +2(t>  tC.n  O+tf 

Q • * Rats  of  mass  flow  lb (M) /sec. 

» . 

R * Radius  in  l y plane  (figure  lit) 

Rt  = Distance  between  point  P and  origin  0 (figure  Hi:) 


Re  = Reynolds  number  - 


V » Absolute  velocity  of  the  fluid  - ft. /sec# 

W * Relative  velocity  of  the  fluid  - ft. /see. 

a * radius  of  the  transformed  circle  in  j ft  plane  (figure  13B) 
b » distance  between  point  on  the  circle  and  origin  0 (figure  13B) 
c * Constant  of  Joukowsky1 s transformation  equation 
i » Imaginary  unit 

1 * component  of  MO  in  J-  direction  in  % of  c (figure  13B) 
m • Distance  between  M and  0 (figure  13B) 
p « pressure  lb(F.)/ft.2 

q*  * k 
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Radial  distance  in  the  cylindrical  coordinates 
Radial  distance  of  the  casing  wall  - ft* 

Radial  distance  of  the  hub"  - ft. 

Radial  distance  of  the  blade  tip  - ft* 
Tip-clearance  - in  (or  ft.) 

v 

Time  - second 

Relative  velocity  of  the  inlet  fluid  - ft/soc* 


Relative  velocity  components  in  r,  y?  and  z direction  respect 
tively  - ft/sec. 


z*  + y* 


- distance  in  the  tangential  direction  - ft* 


Axial  distance  in  the  cylindrical  coordinates 
Horizontal  component  in  yf  - z*  plane  (figure  13/0 


Angle  of  attack  of  the  inlet  relative  fluid  - deg# 
Angle  between  MN  and  £ - axis  (figure  13B)  - deg* 
Angle  between  Mo  and  2 - axis  (figure  13B)  « deg* 

f'i 


* ^ 4-  i>J  (figure  13B) 


Anglo  between  OP  and  J -axis  (figure  1U) 


for  rotor. 


for  stator 


Viscosity  of  fluid  - lb(F)-sec/ft, 

Kinematic  viscosity  of  fluid  * - ft.^/sec* 

Radius  distance  from  toe  blade  tip  inside  the  clearance  space 
- in  (or  ft.) 
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F • 
<r  - 

Cf  - 

y*  - 
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density  of  fluid  - lb,  (l^l) /sec* 

c 


Angular  displacement  in  the  cylindrical  coordinates  - deg. 
Blade  cascado  stagger  angle  - deg. 

Angle  between  PM  and  ^ -axis  (figure  lU)  - deg. 

Rotating  speed  of  the  turbbmachino  - rad,/sec. 


Subscripts 

a ■ At  the  contour  of  the  airfoil  section 
c » Casing  wail 
h ■ Hub 
t c Blade  tip 

r) 

'f  )■  In  the  radial,  tangential  and  axial  directions  of  the  cylindrical 

z) 

0 coordinates  respectively. 


BASIC  EQUATIONS 


The  three-dimensional  flow  of  a viscous  incompressible  fluid 
is  governed  by  the  following  set  of  the  basic  laws  of  fluid  mechanics* 
From  the  principle  of  conservation  of  matter / 

li-  <■  p*-*  a) 


From  Newton* s second  law  of  motion, 

iDy. 

& t 


=-l 7p+VV*V 

i r 


(2) 


or 


zrr  + (v-v)  V « 

i t. 


2 — 


- -j-:v  P + £ V * 


(2a) 


Since  the  absolute  motion  of  fluid  flow  in  a turbomachine  is 
generally  unsteady  while  the  relative  motion  with  respect  to  the 
blade  is  essentially  steady,  the  preceding  equations  are  further 
expressed  in  terms  of  the  relative  velocity  YJ,  which  is  rolatod  to  the 
absolute  velocity  by  the  relation  , 

V =•  /V  + iv  XT  (3) 

In  equation  (3)  l~'  is  the  angular  velocity  of  the  blade  aboyt  the 

J'  -axis)  and  F is  the  radius  vector  measured  from  the  z - axis. 


Now,  for  the  blade  rotating  at  a constant  velocity  o>)  about  the 


z-ajds. 


But 


BB  ~ _ CJ-  p + 2 cO  XW 

o t P t r 


jl  = + ( w-r)  *v  = J~  - M*  (»5{/(W 


(W 


When  the  preceding  relations  arc  used,  equations  (1)  and  (2)  become 
(assuming  steady  relative  flew) , 


— . 
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7 • W ^ o 


(5) 


and 


77W1'  WX^XW  ) -to  f + 2 U>  XW  ~ ~ f &P+P  7 V (6) 


Since  the  boundary  walls  of  a turbomachinc  arc  surfaces  of 
revolution,  a relative  cylindrical  coordinate  system  /",  ^ and  z 
with  measured  with  respect  to  the  rotating  blade  and  rotating  with 
angular  speed  W about  the  z-axis  is  employed.  Then  equation  (5) 
bocomos 

<7> 

and  equation  (6)  gives 

(a)  in  radial  directions  ; 

w,  . f » * - J # - -t*  - « <■  ' 2 ,0  **  , . 

„ UHL  tf(tA  u m * , 2«t  - $ - M^) 

\?t*  *r  cf"  r*  fu*  f zjt*  ^ t a J 


(b)  in  circumferential  direction: 

?iV*  , 2 yjf  . dWu 


(3fl) 


~.±2£  >t?*z 

- er?f  - 


f?*r  ,?wV  I ?Wv 

bf~  VtfY*  ;p  " fr*) 

(8b) 


(c)  in  axial  direction  , 

= -iv<  * W ^ ’t&'rni  'sp 

Cs  V . . « m.  « . /a  \ » 


(8c) 


Although  equations  (7),  (8a),  (8b),  and  (8c)  give  four  indepen- 
dent equations  for  the  "solution  of  four  unknowns,  VV^  vV^  A ^ 

there  is  yet  no  mathematical  methods  available  for  an  exact  solution  of 
theso  oquations.  In  the  following,  the  usual  procedure  of  an  ordetvof- 
magnitude  analysis  is  used  to  obtain  a set  of  simplified  oquations  for 
an  approximate  solution  of  the  problem* 


i 


* *» 


*Y*  *•» 


SIMPLIFIED  EQUATIONS 


We  first  transform  equations  (7)  and  (8)  into  dimensionless 
forms.  Let  tho  representative  velocity  be  W0>  which  is  the  relative 
velocity  of  approach  to  the  blade  just  outside  the  boundary  layer  at 
the  wall.  Let  the  representative  length  be  L,  which  may  be  either  tho 
radial  longth  of  the  blade  or  the  chord  of  the  tip  section  of  the  blade* 
Then  the  equations  arc  made  dimensionless  by  introducing  tho  following 
dimensionless  variables 


Hi  - -*4 


- £- 


r ui  t. 


r Hr 


r 


< - ~r 
* w. 


Equations  (7)  and  (8)  then  become 


r f + 5?'  ° 


r n 


r if  7$ 


(9) 


-iv  r - m/. 


K Pi  f-  T!i.  1«<L  + «:  . »i 

_ A'  i t 

= + ± J ?K'  I }w'  ,?*<.'  -jV//  vv{-'  2 a *5?  ] 

?>'  Keljp:  *pi?  ~ 7* 'P*?  f1**) 


K-p  «- •*  & f < T?r  + JXJ» 


' 

V ■ af 


(V/< 

?' 


'*¥  Kejjfh  }'  ~?r:  *V*  r > *?] 
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#1  + 


|V; 


r 1*f 

= Jii  + ±J  f*4  . ± M , .1  i*& 

W Ke/Tr^  r ?r 


/ 

4 


?H'  l 


where 


K/f^- 


/(10c) 

(11) 


Now,  r»,  z»,<*>  W^>  , Wz*>  p*  arc  of  the  same  order  and  taken 

to  be  of  order  1*  The  dimensionless  tip-clearance  s*  («  s/L)  is 
snail  of  the  first  order,  the  same  as  the  dimensionless  boundary-layer 
thickness  along  the  wall  («*  3/l).  ^hey  belong  to  the  next  smaller 
order  than  1 and  are  indicated  by  0(  £»).  \/\jt  also  has  the  order  of 


Then,  because 


- 


i*r-  Jr  - c(l r j 

' 


»lf V/r/ 


s V'Vz.  / / 

is  of  the  order  1.  On  the  other  hand,  /V^  and 

are  of  the  order  1 (except  at  the  bounding  wall  and 
the  blade  tip  surface) . Based  on  these  orders  of  magnitudes,  we  found 
the  order  of  magnitudes  of  all  other  terms  in  the  quations  and  they  are 
listed  below  j 

X. 

'± 


j-. 


o<jr) 

cl') 


7 , ' 

7>  w'. 


/*• 


Dr-'^ 


and 


c(D 


, 

zr  , 

■a  K// 

# 

' . r 

n* 

.1  ^ 

?w/ 

r’  ?f  > 

H‘ . ' 

' , ?hL 


, ?iy,  1 


/ ?K 


* 
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When  those  orders  of  magnitudes  are  used,  and  a term  of  order 
( &' ) is  omitted,  tho  continuity  equation  (9)  becomes 


Jr7'  f 


— * - o 


(12) 


The  relative  importance  of  the  various  terms  in  equations  (ip) 
depend,  however,  on  the  order  of  magnitude  of  the  Reynolds  number. 


For  fluids  of  small  viscosity,  such  as  air  and  water  at  ordinary 
temperature  and  pressure,  the  characteristic  Reynolds  number  is  usually 
large  and  is  in  the  order  of  (l/jr'^  ).  For  other  fluids,  such  as 
lubricating  oil  and  some  liquid  chemicals,  V is  a hundred  times  greater 
and,  especially  with  relatively  small  characteristic  length  L,  the 
Reynolds  number  is  of  the  order  of  ( ' ) . The  following  simplified 

equations  are  therefore  given,  for  these  two  orders  of  magnitude  of  the 
Reynolds  number  respectively. 

Case  I.  Reynolds  number  of  order  of  \/$ . In  this  case,  equations 
(10)  roduce  to 


/2- 


y£-  = .»i  * IV  rV. 2.  «iW 

?-r  f ' R* 

i j.  z*ik 

r'  ?<f  ~ Re  i}‘* 


2JL  _ JL  V*bL 
?$'  Re 


(13a) 

(13b) 

(13c) 


* i 


In  equation  (13a),  we  see  that  the  viscous  term  is  of  the  same 
order  of  magnitude  as  the  inertia  term.  They  both  contribute  to  the 
radial  prossuro  gradient. 


— - I Rfc. 


r-'  h » 


In  each  of  equations  (13b)  and  (13c),  only  one  viscous  term 
predominates  and  essentially  determines  tho  circumferential  and  axial 
pressure  gradient  respectively. 

It  is  noted  that  the  terms  on  tho  right-hand  side  of  equation 
(13a)  are  of  the  order  of  1,  while  those  on  the  r ight-hand  side  of 
equations  (13b)  and  (13c)  are  of  the  order  of  Therefore,  the 

pressure  gradient  in  the  radial  direction  is  much  smaller  than  the 
pressure  gradients  in  the  othop  two  directions'*  In  other  words,  only 
equations  (13b)  and  (13c)  have  to  be  considered  in  this  case* 


I / * 

Case  II*  Reynolds  number  of  order  of  /<*'  • In  this  case  equation 
(10)  reduces  to 


If  _ (U/ytUJ'r')  _ Vtg* 

1?  ~ r ' >' 

j.ZL  - i 

r iV  ~ Re  " ( ' w 1 t"  n d W} 
2ii  _ j.  ?j&  _ 

%>  Re  ct  r 


(Hia) 

(lUb) 

Uc) 


The  inertia  term  now  predominates  in  the  right-hand  side  of 
equation  (lha),  which  is  recognized  as  the  " simplified-radial-equili- 
brium equation. 


It  is  noted  that  in  this  case  all  terms  that  appear  on  the  right- 
hand  side  of  tho  three  equations  are  of  the  same  order  of  magnitude  of 
1,  Therefore,  the  pressure  gradients  in  the  throe  directions  are  in 
general  of  the  same  ordor  of  I, 
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CASE  I-  WITH: 


-oik) 


Governing  Equations  and  Boundary  Conditions 


With  the  pressure  gradient  in  the  radial  direction  negligible 
compared  to  that  in  the  circumferential  and  axial  directions,  we 
have  the  equations  of  motion  in  the  latter  two  directions  in  dimen- 
sional forms ( 

x 


/ 


u?Jtx.  - 1 lL 

' " r Tf 


?r 


(15a) 

(15b) 


The  simplified  continuity  equation  in  dimensional  form  is 

irt  + ±^  + ?«t 

7r  r ?>fr 

The  boundary  conditions  are  as  follows: 


(16) 


(1)  At  the  blade  tip  surface  i.e.  at  f - , (see  Figs.  8 to  10), 


- Wy,  =:  *o 


(17) 


(2)  At  either  casing  or  hub, 

=-  VV|  =o 

(3)  At  the  casing  for  rotor  blades  (Figs.  8 and  10) 

Wy  » - u)  rc 

At  the  hub  for  stator  blades  (Fig.  9) 

Wy,  rs  lOtj 


(18) 

(1?) 

(19a) 


'W 


*7  **  v i ^ ; * 
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Variation  of  Velocities  in  the  Clearance  Spacing 

Now,  with  pressure  within  the  clearance  space  considered  to  be 
constant  in  the  radial  direction,  the  pressure  p is  a function  of  only 
if  and  z.  So  are  and  . When  equation  (l$a)  is 

integrated  twice  with  respect  to  t at  a fixed  set  of  values  of  ( f 
and  z,  we  obtain  . 

= Jrif  r-(Ur-i)  + c,  ^ +c2 

The  constants  and  C2  are  determined  by  the  use  of  boundary  conditions 
(17),  (19)  and  (19a),  In  order  to  obtain  one  equation  for  either  rotor 
blade  or  stator  blade,  we  let. 


JL 

rt 

Ik 

rt 


for  rotor 
for  stator 


(20a) 

(20b) 


Then. 


fi(A>rt-i)  +c,a  t c 2.  «<? 

+c,Art  = +M*  ft 

where  the  minus  sign  in  the  last  term  is  used  for  the  case  of  rotor 
blades  and  the  plus  sign  is  used  for  the  case  of  stator  blades.  We 
solve  the  preceding  equations  for  C1  and  C2  and  obtain 

*>■  t'-'-r,  ■ ^ ) j, % t-±w 


Hence, 


/>/*>*  /£  \ r 

TIT  (ft.  • > 7 


A-( 


(ZJ) 


* 


«*• 
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Similarly,  integrating  twicef  with  respect  to  f~  , we 


obtain, 


W _ _L  2i 
2 


r + c^r  + c4 


After  the  boundary  conditions  arc  Used  to  determine  the  constants  C- 


3 


and  Cj^,  we  have, 


(22) 


Where  /)-—•>/  for  rotor  blades  and  A «*  Jb.  ^ / for 
stator  blades. 


for  rotor  blades  and  A » / 


Equations  (21)  and  (22)  then  give  the  complete  radial  variation 
of  the  tangential  and  axial  velocities  over  the  clearance  space* 

We  see  that  the  tangential  velocity;  is  composed  of  two  parts.  The 
first  part  is  associated  with  the  circumferential  pressure  gradient  and 
the  second  part  is  associated  with  the  relative  motion  between  the 
blade  and  the  wall.  With  the  coordinate  system  chosen  herein 
(see  Fig, 7 ),  the  circumferential  pressure  gradient  is  a positive 

quantity  either  at  the  rotor  blade  tip  (Fig.  Q ) of  a compressor  or 
at  the  stator  blade  tip  of  the  turbine,  and  is  a negative  quantity 
either  at  the  stator  blade  tip  of  a compressor  (Fig.  ) or  at  the 
rotor  blade  tip  of  a turbine  (Fig.  /o).  The  sign  of  the  multiplier 
of  the  pressure  gradient  can  be  seen  by  the  use  of  series  expansion 
as  follows; 


(!)  Variation  of  Tangential  Velocity  over  the  Rotor  Blade  Tip. 
In  this  case, 


Let 


A = JL 

rt 

L 


fl 


/ V- 

ft 


(23) 


y w >> 
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Then 


A — / -A 
6v ^ =-  ;H(tt 6") 


^ , 3 4 

* 6 " —■  **  , <r* 

* 3 <?. 


(2k) 


(25) 


Let  also 

r 

n*?  i 

4f 

~~  = ' + ~*  = ,+  e 

(26) 

Then 

2.  q?  J-A 

= In  (l  +6  ) * £ J= + ~ + 

Z 3 4 

(27) 

and  equation  (21)  becomes 

it 

6((T-6) 

2 

} , <*jjr  + ^ ^ ^ 

( 3 * ■ i j M ?<f 

Let  ■)) 

s ± 

/ 

S 

(28) 

(Notes  )j 

Z o Of 

r ~ ^ -f  Qf  Aor/x) 

Then, 

i ' ! ' '-  v 

^ = -/1P-fi-  (~1)*  f 


■ ■■)■ 


)>f 


j /<  i 

- >/  (29) 

~ **  1 “*  «•  « ^ than  unity,  th 

qXlty  tnSide  the  Wket  13  ^ * the  coefficient  , 

^ 15  “XWayS  ne8atlTC-  Th«  *«  on  the  right-hand  aide 

of  the  equation  in  eubotantiaUy  proportional  to  the  equare  of  the 

oloaranc,  . and  to  the  circumferentia!  procure  gradient.  It  ie 

identically  equal  to  aero  at  the  blade  tip  and  at  the  eating  eall.  0 

he  other  hand,  the  magnitude  of  the  aecond  ten.  on  the  right-hand  ei 

of  the  equation  inoreaee.  from  zero  at  the  blade  tip  to  *YC  at  the 
casing  wall. 


y .<* 


In  the  case  of  the  rotor  blades  in  a compressor,  the  pressure 
gradient  is  positive  and,  consequently,  both  torihs  on  the  right-hand 
side  are  negative  and  their  magnitudes  add  together. 

In  the  case  of  the  rotor  iblades  in  a turbine,  the  pressure 
gradient  is  negative,  and,  consequently,  the  two  terms  counteract  each 
other. 


The  variation  of  the  coefficient  in  the  first  term  on  the  right 
hand  of  equation  (21)  is  plotted  for  several  values  of  $ in  Fig,  // 
It  is  seen  that  its  magnitude  increases  at  first  with  radius,  reaches 
a maximum  at  equals  to  about  o,s  and  then  decreases  to  zero 
at  the  casing. 


(2)  Variation  of  tagential  velocity  over  the  stator  blade  tip. 
In  this  case. 


x - 4 = — * /-«• 


/i 


In  a = -<*+§ 


_/-  _ BjiJ  ~ i~  & 


rt 

r 


n 


y3-  i .4 
6 b £ 


In  — - - (*  + z_  + f + 4 


(21*) 

(25a) 

(26a) 

(27a) 


and  equation  (21)  becomes 

//0V ) } 6 +<j'f'  tt  7 ft Z'l  f-  ~ 

r-+—r*"  17~?r  r 


T 


- 2 ' ->  r j r 


/ 3.  <n  J ^ ~*J 

We  Bee  again  that  the  quantity  inside  the  bracket  is  always  a 


positive  one  and  that  it  is  always  slightly  greater  than  that  in 


equation  (29)  for  the  same  value  of  <$  and  >J  • In  the  case  of 
a compressor,  the  pressure  gradient  now  is  negative,  both  terms  on  the 
right-hand  side  of  equation  (30)  are  positive.  In  the  case  of  a 
turbine,  the  pressure  gradient  id  positive,  and  consequently  the  terms 
counteract  each  other.  The  first  term  varies  with  y in  a similar 
manner  as  that  in  the  rotor  case.  The  second  term  varies  linearly 
from  a value  of  zero  at  rt  to  a'  value  of  tofy  at  • 

The  variation  of  the  axial’  velocity  as  given  by  equation  (22)  can 
also  be  put  in  terms  of  the  dimentionless  variable  J for  either  the 
rotor  or  the  stator  blade  as  follows: 

W,  = - yU-},)j£*L  (31) 

Sr  t ( 7$ 

We  see  that  the  axial  velocity  is . proportional  to  -sz  and  the  pressure 
gradient  and  it  has  a parabolic  variation  with  ^ inside  the  clearance 
space.  Since  the  pressure  gradient  in  the  axial  direction  is  always 
positive  for  compressor,  the  axial  velocity  is  always  negative; in  the 
case  of  a turbine,  the  axial  pressure  gradient  is  negative  and  the  axial 
velocity  is  positive. 


Variation  of  Pressure  in  the  Clearance  Space 


In  order  to  compute  the  velocity  variation  by  equations  (21)  and 
(22)  throughout  the  clearance  space,  it  is  necessary  to  determine  first 
the  pressure  gradient  throughout  the  space.  This  is  done  as  follows: 


Multiplying  the  continuity  equation  by  e/f~  and  integrating 
from  the  tip  of  a rotor  blade  to  the  casing  gives , 

f w dr  t r%£  + I dr  * o 

J,-t  * k ' 'rt 


y •%» 


* «*  «!»• 
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But 
Also 
Now,  let 
Then 

Since 


r rc  !r<- 

f dr  = w,. 

\ ” u 

f r,i  j,  = 0*  ^ 


k ' >* 

r'  * /nT 

r U 


k T* 


f.  ej  ur 


( ft  , ^ u/  2>rt  „ / v«e 

| ***  d In  = j J&.jf  - *Y<t*  + tyrtW 

ft  ?if 


X 


Therefore 


iVy;  SO  2J5 


?£■ 

/J'Jf 


■« 


5=  O 


Similarly, 


[ up dr  ,&/  «¥*r  = 
Jrt  ^ ’P 


**  dl- 
l' 


L 


H 


Since 


L^i  , 

v;  ** 


'0  n't  id  r —- 

v~j,dr-  t Wj,r*Jy  *n  if 


Hence 


r f~C 


~ . t-t  ~° 

rc 

i 

The  integral  form  of  the  continuity  equation  is  then 


'»• t 


dr  -£-Jr  %'  e/r 


f m. 

hi 


r< 


»L  ~F  d‘  +?Ht  *>*** 


(32) 


When  the  values  of  |Y^.  and  <V^,  as  given  by  equations  (21) 

A ®”bstituted  into  equation  (32),  we  have, 

-,)_j t/vr]tii}.  , A-ff-nS-ldr 
/ ii  vnhahf  *-<  '*  r I 


V 


■ v 

/ 


MV 


?v  „ 


Ts  **  ■■  >,  (*> 


Carrying  out  integration  and  differentiation,  and  Rearranging,  we 


finally  obtain, 


*(***?)  2?  ?J> 

(33) 


This  is  a linear  second-order  partial  differential  equation  and  can 
easily  be  transformed  into  a Laplace  equation  by  ]e  tting 

u = t rt¥  <*> 

“ j 12  ((*-!)-  M/hA)1-]  / * 

Thus,  j 

Kt 

>Y  * ° <35) 

The  relation  between  ^ and  p can  be  more  clearly  seen  by 
substituting  equations  (23)  to  (27)  into  equation  (3U)*  This  results 

ta  jf  A. i ^ — inr 

v < /-  <f  + 4, to  or  / - - - J 


^ <r  * - • )rt<f  ^ r<f 


(3Ua) 


Thus,  y has  the  same  sign  as  <f  and  is  approximately  equal  to  the  casing 
radius  times  y> 

If  the  continuity  equation  is  integrated  from  the  stator  blade 
tip  to  the  hub,  we  again  obtain  equations  (33)  to  (35)  • But,  when 
the  relations  (23a)  to  (27a)  are  used,  we  obtain  / 
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U =/ L r lrtlf 

o ( /ts  f 4jf  a + - - • • J 

= (/-  <r ;fcf 

~ br 

A general  approximate  form  for  equation  (3i|)  is 


(31 i) 


(36) 


and  an  approximate  form  of  equation  (33)  is 

*■  2 

-L-  JJL  , Zi_  i 0 

/l*  ifVy1- 


(37) 


Equation  (37)  is  a second-order  partial  differential  equation 
of  the  elliptic  type.  Hence,  the  proper  boundary  condition  is  the 
pressure  distribution  around  the  blade  surface  at  the  blade  tip* 

Like  other  boundary-layer  problems,  the  potential  flow  outside 

the  boundary  layer  has  to  be  known  in  order  to  solve  the  boundary-layer 

flow. 

The  pressure  distribution  around  the  blade  at  the  tip  section 
is  mainly  determined  by  the  geometry  of  the  blade  (the  shape,  the 
stagger  angle,  and  the  solidity),  the  angle  of  attack,  the  Mach  number, 
and  the  Reynolds  number.  The  influence  of  these  factors  on  the  tip- 
clearance  flow  is  therefore  exerted  through  the  boundary  condition  for 
the  solution  of  equation  (33)  or  (35)  • As  the  first  approximate 
solution,  the  effect  of  these  factors  on  the  pressure  distribution 
around  the  blade  tip  section  can  be  obtained  either  by  using  the 
theorectical  calculation  based  on  two-dimensional  potential  flow 
around  the  blade  tip  section  or  by  using  cascade  tests.  In  the  actual 
three-dimensional  flow,  the  pressure  distribution  around  the  blade 
tip  section  is  further  influenced  by  the  ratio  of  clearance  to  blade 


I 


height,  the  relative  speed  between  the  blade  and  the  wall,  the  secon- 
dary boundary-layer  flow,  across  the  channel  formed  by  the  blades,  and 

1 

the  three-dimensional  geometrical  shape  of  the  blade.  These  influences 
can  only  bo  determined,  at  the  present  time. by  experiments,  some  of 
which  are  scheduled  in  the  secaind  phase  of  this  research  project. 

With  a given  geometry  and  the  pressure  distribution  around  the 
blade  at  the  tip  section,  equation  (35)  can  be  solved,  and  the  circum- 
ferential and  axial  derivatives  used  in  the  preceding  equations  for  the 
computation  of  the  velocity  components.  The  method  of  solution  will 
be  discussed  after  similar  equations  like  those  contained  in  this 
section  are  obtained  for  the  second  case. 
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CASE  II  WITH  A?tf  * ofW) 


In  this  case  the  inertia  forces  are  in  general  of  the  same 
order  of  magnitude  as  the  viscous  force,  and  the  pressure  gradients  in 

all  three  directions  are  of  the  saicje  order  of  1.  In  dimensional  forms, 

/ 

the  three  simplified  equations  of  potion  are: 


(Vy 'if+U  t-f  _ V * 1 

f r ^ ...  ■* 

ivr  + —£  -p  + = - r rt<f  *r' 


> 


IV  224  + l*&  =-i-\ f 

Tr  +r2<p  p 21  f ?A  7>t 


(38a) 


(38b) 


(38c) 


The  simplified  continuity  equation  is  the  same  as  in  the  previous 


?/  y /■  n ■ if  ~ 


Although  the  radial  pressure  gradient  is  now  of  the  same  order 
of  magnitude  of  the  pressure  gradient  in  the  other  two  directions,  the 
total  variation  of  pressure  across  the  clearance  space  of  ordinary  size 
is  still  relatively  small  (in  the  order  of  two  to  three  percent),  (in 
reference  20,  a measured  difference  of  % across  the  boundary  layer 
along  the  wall  of  a compressor  is  reported).  Therefore,  for  an  approxi- 
mate solution,  we  may  still  ignore  the  radial  variation  of  pressure  for 
our  problem,  i.e.  the  pressure  is  considered  to  be  a function  of  only 
and  z. 

Equations  (38b)  and(38c)  each  have  three  more  inertia  terms  on 
the  right-hand  side  of  the  equation  than  the  corresponding  equations 
of  the  previous  case.  At  either  thq  blade  tip  or  at  the  bounding  walls, 
however,  these  inertia  terms  are  equal  to  zero,  because  either  the 


*•  \ 
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velocity  component  or  its  circumferential  derivative  is  equal  to  zero* 
These  terms  exert  their  influence- at  in-between  radii.  Because  of 
their  presence , the  equations  become  non-linear  and  it  is  found  that 
an  exact  solution  of  equations  (38b),  (38c),  and  (16)  cannot  be 
obtained.  In  the  following  an  approximate  solution  is  obtained  by 
assuming  polynomial  variations  for  the  velocity  components  and 
determining  the  constants  by  the  boundary  conditions  on  the  velocities 
and  their  radial  derivatives  as  given  by  equations  (38b)  and  (38c). 


Thus  we  assume 


c Ac  + A,  I!  f-  A z >j  t A 5 ^ 4-  - - - - (39) 
- Be  +B,yt  y i-Btf*---  m 


The  boundary  conditions  are  as  follows: 
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(lilb) 
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at  ^ / ( r ■ rc  for  rotor  blades  or  A * /J  for 


stator  blades) 


I *Jr  * o 


---  i-  A IV  It 


(ad) 


(ae) 


(where  the  minus  and  positive  signs  are  used  for  rotor  and  stator 


respectively) 


(V  ? - O 


(ULC) 


at  "h  x o , Equation  (38b)  gives 


y ••  ■■ 


Then , 


"V  * 
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??!*- 


-32- 

_i.  iL 

/Mt 

_sl  it 

/'/*  ?f 


(Uig) 


Also,  Equation  (38c)  gives 


/ */?. 


1 it 
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* = ? 


or 
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At  >j  ■ 1,  equations  (38b)  and  (39c)  give,  respectively 
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~''7nyZ  ~ Mfh  ? k 

and  -»*W>  Ipf 
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(ah) 

(ai) 

(103) 
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From  equations  (Ulh)  and  (aj)  we  see  that  the  second  partial 
derivative  of  Wz  with  respect  to  V at  the  blade  tip  ( >7  ■ 0)  is 
the  same  as  that  at  the  wall  ( >7»  1) , In  othor  words,  the  second 
partial  derivative  should  not  include  y , Thus,  equation  (UO)  takes 
the  form , 

= 3C  + t e*  *jz 

By  the  use  of  equations  (ac),  (af),  and  (ah),  we  obtain 

B0  *0 

vr  iL 
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Then 


W„ 
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On  the  other  hand,  equations  (ag)  and  (ai)  show  that  the 


second  partial  derivative  of  W y>  with  respect  to  at  the  blade  tip 
( >j  • 0)  is  slightly  different  from  that  at  the  wall  (^  ■ l)i  Sin6e 
we  have  altogether  four  boundary  conditions  for  VI  if  and  its  derivative 
to  satisfy  at  the  two  end  points,  we  take, 

'J 

W*  = A,  (ii3) 


From  equations  (i*lb),  (l*le),  (l*lg),  and  (Itlj),  we  find 

Ac 

A,  t a,  + A*  K r Aa)r* 
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_s~  it 
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45 
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Solving  for  A0,  A^,  A 2 , and  A#  and  substituting  them  into  (1*3) 


gives, 


— 


'^7 


When  equations  (1*2)  and  (UU)  obtained  herein  are  compared  to 
respectively,  equations  (31)  and  (21)  of  the  previous  case,  we  see 
immediately  that  the  axial  velocities  as  given  by  equations  (1*2)  and 
(31)  are  exactly  the  same.  To  compare  the  tangential  velocities,  we 
substitute  equations  (2i*)  and  (2l*a)  into  equation  (1*1*)  resulting  in 


for  stator  blades.  In  these  forms,  we  see  that  equations  (1*5)  and  (1*6) 
are  also  equal  to  the  corresponding  equations  (29)  and  (30)  of  the 

previous  case  if  terms  involving  c>*  and  higher  orders  are  neglected* 


-3k 

Thie  is  usually  permissible,  since  < is  a very  small  quantity. 

Hence  we  see  that  if  the  radial  variations  of  the  velocity 
components  in  this  case  can  be  represented  by  polynomials,  the 
equations  for  the  velocities  are  the  same  as  those  of  the  previous 
case.  The  validity  of  the  assumption  of  polynomial  representation 
in  this  case  should  be  verified  in  the  future  experiments.  It  should 

i>  | 

be  noted  that,  although  the  equations  are  the  same,  the  pressure  dis-  1 

tribution  around  the  tip  blade  section  and,  consequently,  the  circum-  I 

ferential  and  axial  gradients  of  pressure  can  be  a function  of  Reynolds 
number.  Therefore,  the  actual  velocities  can  still  be  a function  of 
Reynolds  numbers  The  values  of  in  the  two  cases  are  also  different* 

Inasmuch  as  the  velocities  of  this  case  are  the  same  as  those 
in  the  previous  case,  equations  (31*) » (35) > (36),  and  (37)  obtained  in  j 

the  previous  case  can  also  be  used  for  this  case,  I 


in 
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SOLUTION  OF  PRESSURE  DISTRIBUTION  IN 
CLEARANCE  SPACE 


From  the  preceding  two  sections  we  see  that  the  pressure  distri- 

i * 

bution  across  the  clearance  spabe  is  to  be  determined  by  the  equation 

i w 


or 
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(37) 
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(35) 


where 
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for  rotor  blades 
for  stator  blade 


(36) 

(31*0 

(3lib) 


Now,  the  difficulty  in  solving  the  pressure  distribution  lies  in 
the  arbitrary  boundary  values  of  the  pressure  given  around  the  blade 
section  of  agiven  arbitrary  shape*  In  general,  a sufficiently  accurate 
answer  can  be  obtained  by  the  use  of  relaxation  or  matrix  method 
(Reference  20,  21),  The  difficulty  introduced  by  the  non-uniform  grid 
spacings  can  be  helped  out  by  the  use  of  differentiation  formulas  for 
non-uniform  spacings  (Reference  20.  See  example  on  pp.  27-29  for  the 
use  of  these  formulas  in  conjunction  with  either  the  relaxation  or 
matrix  method), 

A second  method  of  solution  is  to  replace  the  given  blade  section  by 
a Joukowsky  airfoil  section  and  to  solve  for  the  pressure  variation  over 
the  airfoil  section  after  a conformal  transformation  which  converts  the 
airfoil  into  a circle.  This  method  of  solution  is  given  below. 


Approximation  of ; the  Given  Tip  Section 
byn  Joukowsky  Airfoil 


Figure  13  shows  the  blade  tip  section  in  the  y - z plane*  JL  is  the 
angle  between  z-axis  and  the  chord  line  which  connects  the  leading  edge 
and  the  trailing  edge  of  the  airfoil.  Also  shown  is  the  Joukowsky 
airfoil  approximating  this  tip  section*  This  Joukowsky  airfoil  is 
obtained  by  matching  the  chord,  ,;the  thickness  ratio  and  the  camber 
(maximum  value  of  mean  line  ordinate). 

Referring  to  Figure  13,  we  first  make  a transformation  of  coordin- 
ates from  y - z to  y«  - z*  by, 

l.  (H7) 

'I'  * J '*"/  J 


Then  this  Joukowsky  airfoil 'with  its  chord  (C)  lying  on  the  z>*  -axis 
is  obtained  by  transforming  a circle  of  radius  >a>  with  center  at  M in 
the  2 - }j  plane  by  the  Joukowsky  transformation  : 

* = ? + j1  m 

where 


*•  * 


(h?) 

(5o) 


c =4c('  + ™z) 

CL  a C {/  + ">) 


(51) 

(5a) 


(m  is  a small  quantity) 
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The  coordinates  of  the  airfoil  are  given  by 


z.  * - ( b + ~ - 2 c-  i ose  ') 

a K b / fiV; 

jj'  = ( b - /t*‘ n 6 * 2.  c *^,6 [M(i+u&e)+p>***>6y 


where 


±‘  n)  *<&*  ~ 

-tnn  — - — — 7^  •v 

c ■*  r>i  sin  « 


rr)  cos  TT 
c.  ■+ msmf 


and  is  a small  angle.  ■ 

From  eqs.  (53) » the  maximum  thickness  is  equal  to  -3  ,7)  C 

which  occurs  at  one-quarter  of  the  chord  from  the  leading  edge  and  the 
camber  is  equal  to  1 /Z  • In  the  approximation  of  the  given  blade 

tip  section  by  the  Joukowsky  airfoil,  the  geometrical  data  of  the  tip 
section  is  to  be  used  to  obtain  C (or  c),  m,  and  ft  of  the  Joukowsky 


airfoil. 


Velocity  and  Pressure  Distributions  Around  the 
Joukowsky  Airfoil 


In  the  ^ -plane,  the  velocity  at  the  surface  of  the  circle  is 


/f  = c2  ( o{  + : 

7 ' ' 2T 


27TC\ 


p - 4-irQ  f/3) 

Hence 

* 3 Wo  * £)  + ^ *i(o(+p )J  ( 5'5* ) 

The  velocity  around  the  airfoil  in  the  / -plane  and  that  around 
the  circle  in  the  kj  -plane  are  related  by 
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where 
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From  Fig.  13,  we  see 


b sh'-yi  6 — Q.  (/oU\  -+-  $ ,.‘V1  ^ 

tcs  0 **  cl  ocs  Tjr  +■  / c 


(58) 
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Hence 
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Substituting  equations  (55),  (57)  and  (60)  into  equation  (56)  gives 

f ^ 
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Equation  (6l)  now  gives  the  velocity  distribution  around  the  airfoil 


in  terms  of  \l 


o . 


o< 


p,  * 


and  . For  a 


given  airfoil  and  a given  approach  velocity  at  a given  angle  of  attack, 
ft.  ft , *V0  and  oi  arc  fixed  values, 

\ 

The  pressure  distrioution  around  the  airfoil  is  obtained  by  the  use 
of  Bernoulli  equation  and  the  velocity  distribution  given  by  Equation 
(61) 


P *' 
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Solution  of  the  Laplace  Equation  in  Pressure 


With  the  boundary  values  of  the  pressuro  around  the  airfoil  section 
given  by  equation  (62),  the  variation  of  pressure  across  the  blade  tip 
is  obtained  by  the  solution  of  equation  (35) • Because  of  the  complica- 
ted boundary  shape  of  the  airfoil  in  the  plane,  we  attempt  to 

solve  the  equation  in  conformall- transformed  < -plane  with  the  simple 
circular  boundary.  To  do  this,  the  Laplace  equation  is  expressed  in 
polar  coordinates  R and  P (Fig*  lit)  • 

JJl  .12L  + 1±?-  =0 

iflz  /?  2/?  /?*  (63) 


The  solution  of  this  equation  can  be  expressed  in  a Fourier  series 
as  follows  (see  reference  22) 

f(&,  # ) ~ A0  i-  (6W 


In  order  to  determine  the  Fourier  coefficients  A*s  and  Bfs,  we  first 
write  equation  (6U)  at  the  boundary  R * as 

^ ) = Ac  ♦ 0 (An  + Bh  (65) 


Then  wo  express  tho  known  boundary  values,  equation  (62)  in  a Fourier 
series: 


p (a,  ty  ) - a'  + Jr  ( oasny  f B*n  AXn7if) 
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Equating  equations  (65)  and  (66)  gives 


Ac  - Ao 

An  - A„ 


\ 


(67) 
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In  the  case  of  cambered  Joukowsky  airfoil,  it  is  quite  a complicated 
process  to  express  explicitly  the  A 5 and  B 's  in  terms  of  W 

C 

/ 

(5  and  jf  . 


In  case  it  is  desirable  to  solve  the  present  problem  with  an 
experimentally  measured  pressuro  distribution  around  the  tip  blade 
section,  one  can  obtain  the  coefficients  $ and  B„'.5  by  a 
suitable  numerical  method  and  then  obtain  A n s and  by 

equations  (67)* 


In  the  case  of  a symmetrical  Joukowsky  section,  it  is  found 
relatively  simple  to  obtain  a general  expression  for  the  coefficients 
(based  on  the  theoretical  pressure  distribution.  In  this  special  case. 
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and 
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Based  on  this  equation  we  found,  by  the  method  given  in  reference  23 
and  using  30°  intervals, 
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These  expressions  are  general  that  it  can  bo  used  for  any  values  of 
^ and  c*  • 


Pressure  Gradient 


To  obtain  the  pressure  gradient,  we  differentiate  equation  (6ii)  with 


respect  to  R and  ^ , respectively 
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To  lead  these  derivatives  to  and  JT 
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, we  have  to  find  the 


Prom  equa- 


relations  among  ^ .21  J W •# 

tions(53) j (5h), (55)  and  referring %o  Figfure  13,  *we  obtain, 
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Equations  (71)  and  (72)  give, 
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Differentiating  the  preceding  equations  with  respect  to  z*  gives, 
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By  sol  ing  for 
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we  obtain 
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Similarly,  differentiating  with  respect  to  y»  leads  to 
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In  particular  at  R * a ■ ( / +/ ) c 4*<. ^ , equations  (73)  and  (7U) 
reduce  to, 
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Where,  ;i  * , . 
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Having  found  those  derivative^,  we  can  compute  / ) ane^ 
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Where  and  ( can  be  from  equation  (70)  by- 

putting  R ■ a, 
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From  Equation  (36) 


Equations  (70),  (73)  > (?U),  (76)*"  (77)  and  (78)  altogether  give 

, which  are  reqjxLrod  to  confute  the  velocity  distri- 

v 

bution  across  the  tip  clearance  space  as  given  by  equations  (1|2)  and 

m or  (U6)  .*  ' ; 

Alternate  Approximation 

Instead  of  approximating  the  given  tip  blade  section  by  a 
Joukowsky  section,  as  stated  oh  p/36,  one  may  transform  the  given 
blade  section  by  the  Joukowsky  transformation  into  a near  circle, 
which  can  then  be  approximate  by  a true  circle,  thereby  determining 
the  values  of  rrj  and  /3  • The  rest  of  the  procedure  remains  the 
same. 


and 
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mss  FLOW  ACROSS  TIP-CLEARANCE  SPACE 


After  the  velocity  distribution  across  the  tip-clearance 
space  is  obtained,  the  mass  flow  across  it  can  be  obtained  by  adding  ' 
together  the  mass  flow  in  the  tangential  and  in  the  axial  directions. 
In  the  tangential  direction,  the  mass  flow  is  (Figure  1$) • • 


■ t O c A 
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In  the  axial  direction,  the  mass  flow  is, 
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Therefore,  the  resultant  mass  flow  across  the  tip-clearance  space 
will  be,  (Figure  1$) 
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A basic  analysis  is  made  on  the  flow  phenomena  across  the 
blade  tip  clearance  space.'  The  fluid  is  considered  to  be  viscous 
and  incompressible  and  the  basic  equations  of  motion  and  continuity 
are  expressed  with  respect  to  a cylindrical  coordinate  system  rota- 
ting with  the  blades  at  the  same  speed; 

Simplified  equations  are  then  obtained  by  an  order-of-magnitude 
analysis.  Two  cases  are  considered;  one  with  low  Reynolds  number  and 
one  with  high  Reynolds  number. 

In  the  first  case,  the  equations  of  motion  can  be  solved  to 
give  an  analytical  expression  for  the  velocity  distributions  in  the 
clearance  space  in  terms  of  the  pressure  gradient  and  the  relative 
speed  between  the  blade  and  the  bounding  wall*  The  pressure  distri- 
bution is  governed  by  a Laplace  equation  which  is  obtained  from  the 
continuity  equation* 

In  the  second  case,  the  equations  of  motion  cannot  be  solved 
analytically.  An  approximate  solution  is  obtained  by  assuming  a third 
order  polynomial  for  the  tangential  velocity  and  a second  order  poly- 
nomial for  the  axial  velocity  and  using  the  proper  boundary  conditions* 
The  form  of  the  equations  obtainod  is  very  similar  to  that  of  the  first 
case* 

Three  methods  of  solution  are  discussed  for  the  solution  of  the 
Laplace  equation  of  prossure  distribution.  One  is  a direct  numerical 
solution  by  either  the  relaxation  method  or  matrix  methsd.  The  next 
one  involves  an  approximation  of  the  tip  section  by  a Joukowsky  section 
Tho  last  one  involvos  a transformation  of  the  given  blade  section  by 


the  Joukowsky  transformation  into  a near  circle,  which  is  then 
approximated  by  a true  circle*  In  either  of  the  last  two  methods  the 
pressure  distribution  in  the  transformed  plane  is  expressed  in  the  form 
of  a Fourier  series.  Either  the  theoretical  or  experimentally  obtained 
pressure  distribution  along  the  blade  surface  can  be  used  as  the 
boundary  values.  Expressions  are  obtained  to  convert  the  pressure 
gradients  obtained  in  the  transformed  plane  to  those  in  the  actual  blade 
tip  section. 

Finally  an  expression  for  the  mass  flow  across  the  tip  clearance 
space  is  obtained* 

The  pressure  gradient  over  the  blade  section  at  the  blade  tip  is 
found  to  be  affected  by  factors  such  as  blade  loading  (i,e».  pressure 
distribution  over  blade  section)  shape  of  blade  airfoil  section, 
casing  wall  boundary  layer,  viscosity  of  the  fluid,  relative  motion 
between  the  casing  wall  and  the  blade  tip,,  Velocities  through  the  tip- 
clearance  is  found  to  be  primarily  determined  by  the  pressure  gradient 
and  varies  with  the  square  of  the  tip-clearance,  while  mass  flow  across 
the  tip-clearance  is  found  to  be  determined  by  the  velocity  distribution, 
hence  by  the  pressure  gradient  and  varies  with  the  cube  of  the  tip- 
clearance 
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FIG.  1 - Flow  Past  Free  Tip  Blade. 


no.  2 - Flow  Past  Blades  with  Bad*  Wall. 
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FIG.  3 - Flow  Fast  Bladoa  With  Stationary  End- Wall  And 
Tip-Clearance. 
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Flow  Fast  Blades  With  Moving  End-Wall  And 
Tip-Clearance  in  a Compressor. 
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FI (3.  4A  - Details  of  Flow  Inside  Clearance  Space 


Flow  Past  Blades  With  Koving  Bad- trail  And 
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13/v  ~ Joukcveky  Airfoil  in  F3G„  13B  * Joukowsky  Airfoil  Transformed  into 
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